Tunable terahertz emission from difference-frequency in biased superlattices 
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The terahertz emission from difference-frequency in biased superlattices is calculated with the 
excitonic effect included. Owing to the doubly resonant condition and the excitonic enhancement, 
the typical susceptibility is larger than 10 -5 m/V. The doubly resonant condition can always be 
realized by adjusting the bias voltage and the laser frequencies, thus the in-situ tunable emission 
is efficient in the range of 0.5-6 terahertz. Continuous wave operation with 1% quantum efficiency 
and fiW output power is feasible while the signal absorption in undoped superlattices is negligible. 

PACS numbers: 78.20.Bh, 42.65.An, 42.65.Yj, 73.21.Cd 



Terahertz (THz) electromagnetic waves have poten- 
tial applications in many fields, like medical diagnosis, 
environment monitoring, high-speed communication, as- 
tronomy spectroscopy, etc. Though intense tun- 
able THz emission from free-electron lasers is avail- 
able in laboratories Q, tunable tabletop THz sources 
are still desirable for practical applications. To obtain 
THz emission from small semiconductor devices, many 
schemes have been studied, such as coherent phonons 3| , 
wave-packet oscillation in asymmetric quantum wells ^ , 
heavy-light hole beatings [fj, Bloch oscillations Q, and 
difference- frequency in doped quantum wells Q. With 
state-of-the-art design of superlattices, a prototype of 
quantum-cascade THz lasers has been demonstrated re- 
cently Among these mechanisms for THz emission, 
the difference-frequency process is of special interest be- 
cause of its in-situ tunability, intense output under phase- 
matching condition, and flexibility of operating at both 
continuous wave and pulse modes [9j . Furthermore, dou- 
bly resonant condition, in which both the input and out- 
put are near resonant with transitions in the system, can 
also be exploited to enhance the difference-frequency [j| • 

Doubly resonant difference-frequency in biased super- 
lattices was also proposed for THz emission |lOl 111) . Un- 
der doubly resonant condition, the Bloch oscillation is 
sustained and amplified by the effective THz potential 
resulting from the dipole interaction of excitons and the 
bichromatic input light, generating efficient THz radia- 
tion Several advantages of this method over other 
difference-frequency schemes can be expected: First, the 
applied electric field breaks the inversion symmetry of 
the system, leading to a large intraband dipole matrix 
element. Secondly, the doubly resonant condition can 
always be accomplished by adjusting the static electric 
field and tuning the input light. And thirdly, the prob- 
lem of signal absorption can also be avoided in undoped 
superlattices. 

Though it has been well-known that the exciton corre- 
lation plays an essential role in THz emission from Bloch 
oscillation in optically excited superlattices [12j| . its ef- 
fect on difference-frequency in biased superlattices is still 
unclear. This question will be addressed in this Letter, 



and it will be shown that the excitonic effect can enhance 
the emission power by at least two orders of magnitude, 
which, however, is absent in, e.g., difference- frequency in 
doped quantum wells [7j . 

The second-order difference-frequency susceptibility is 
the key quantity determining the emission intensity. In 
principle, it can be evaluated from the textbook for- 
mula derived with the double-line Feynman diagrams 
fl3| . Under the doubly resonant condition, the difference- 
frequency susceptibility [LJ 
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where fli is the frequency of the input light polarized 
at ej i direction, e aa > (a, a' = a, b, or 0) is the transi- 
tion energy between the exciton states |a), \b), and the 
semiconductor ground state |0), d Qa / is the dipole ma- 
trix element, 72 and 71 are the interband and intraband 
dephasing rates, respectively, V is the volume of the sam- 
ple, and eo is the vacuum dielectric constant. 

With the excitonic effect neglected, the susceptibility 
of biased superlattices can be analytically evaluated , 
and the result turns out comparable to that of doped 
quantum wells and larger by many orders of magnitude 
than that of bulk semiconductors [13L ITij . The Coulomb 
coupling, however, makes it a formidable task to calcu- 
late the susceptibility directly from Eq. , since all the 
excitonic eigen states should be obtained. To avoid such 
an exhausting work, we have developed a time-domain 
technique, in which the susceptibility is first transformed 
into the time domain, and the result numerically calcu- 
lated is transformed back to the frequency domain by 
standard fast Fourier transformation. 

By Fourier transformation of Eq. , the time-domain 
susceptibility of biased superlattices can be derived as 
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in which the exciton Green's function satisfies the motion 
equation in the tight-binding model 

ind t ^( P ,i,t) = [E g - iPj(p,i,t) 

+ J2 {leFD5 mfi - A m /4) ^(p, I + m, t) 

m 
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where X rn is the tunnelling coefficient between quantum 
wells separated by m barriers, p denotes the in-plane 
coordinates in real space, I is the index of the unit cell 
of the superlattice, E g is the distance between the cen- 
ters of the electron and hole minibands, /i is the reduced 
effective mass of the electron-hole pair, D is the super- 
lattice period, F is the strength of the static field, and 
V(p, I) is the Coulomb potential. In the derivation of Eq. 
(J5J), only the lowest electron and heavy-hole minibands 
are included, and the Coulomb potential, assumed slow- 
varying as compared to the superlattice potential, takes 
the form 

V(p,l) = -e 2 (4 7 re 6)- 1 (p 2 + l 2 D 2 y 1/2 , (4) 

where e is the dielectric constant of the material. It would 
not be difficult to include more complexity of realistic 
semiconductor systems, such as the valence-band-mixing 
and the Coulomb coupling between minibands, which, 
however, is expected to modify the results only in details. 

The most important feature of Eq. is that the 
susceptibility has the form of exciton-exciton correla- 
tion. Thus the difference-frequency susceptibility, as well 
as linear absorption spectra, can be evaluated by just 
integrating the motion equation for the exciton wave- 
function [Eq. ©], which can be done with the space- 
time difference method proposed by Glutsch et at (lif . 
To check the numerical method, the susceptibility has 
been numerically calculated with the Coulomb potential 
artificially switched off and compared to the analytical 
result [llj, the deviation is always less than one percent. 

Now let us focus on the excitonic case. To be spe- 
cific, here we consider a superlattice sample that has 
been well studied for Bloch oscillation |l6(, namely, a 
GaAs/Al .3Gao.7As superlattice with 67 A well width 
and 17 A barrier width. A Kronig-Penney calculation 
shows that the lowest electron and hole minibands have 
almost perfect cosinusoid dispersions and the combined 
miniband width is about 41 meV, so A TO rs Si m \ i x 41 
meV. With the miniband width larger than the emis- 
sion threshold of optical phonons, the relaxation in this 
sample is very rapid, so the intcrband and intraband de- 
phasing rates are chosen quite large values: 7-f 1 = 0.6 ps 
and 7^~ 1 = 0.3 ps. The interband dipole matrix element 
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FIG. 1: (a) The 3D plot of the difference-frequency suscep- 
tibility of a biased superlattice vs. the output frequency 
and one of the input frequencies, (b) The linear absorption 
spectrum of the superlattice. (c) The contour plot of the 
difference- frequency susceptibility. 

obtained from the k • p theory is d cv s» 6.5 eA. Other 
parameters are such that the excitonic binding energy is 
4.9 meV, the band gap E g = 1.511 eV, the static dielec- 
tric constant e = 12.9, and the optical refractive index 
n = 3.26. 

The difference-frequency susceptibility has been calcu- 
lated for various static field strength. A typical example 
is plotted in Fig. ^ f° r F ~ H-9 kV/cm (correspond- 
ingly, the frequency of the free-particle Bloch oscillation 
i^o = 2.44 THz, or hv BO = eFD = 10 meV). The dou- 
bly resonant effect is evident in the peak features of the 
susceptibility spectrum. The Coulomb interaction renor- 
malizes the interband transition energy and the Bloch 
oscillation frequency. More importantly, as compared to 
the free-particle result [see Fig. |U(a)], the excitonic ef- 
fect enhances the susceptibility by more than one order 
of magnitude. The excitonic enhancement results mainly 
from the enhancement of the oscillator strength at band 
edge (due to the Sommerfield factor). The Sommerfield 
factor also induces extra absorption at band edge, but 
considering the fact that the linear optical absorption is 
proportional to the second power of the interband dipole 
matrix element while the THz signal intensity is propor- 
tional to the fourth power, the net effect of Coulomb 
interaction is advantageous to the difference-frequency 
process. 

To estimate the power of the THz emission, we assume 
perfect phase-matching for the optical mixing, which can 
be achieved just geometrically since the refractive index 
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FIG. 2: (a) The difference-frequency susceptibility and (b) 
the power of the THz emission as functions of the frequency 
difference. In (b), the THz field strength F^c is also indicated. 
The solid lines with open circles represent excitonic results at 
fixed field strength F = 11.9 kV/cm and optimized input light 
frequency, while the solid lines with close circles correspond 
to excitonic results with both the static field and the optical 
frequencies adjusted to maximize the emission power. The 
dotted lines plot the free-particle results at fixed electric field 
(11.9 kV/cm) and optimized optical frequencies. For visibil- 
ity, the susceptibility and emission power in the free-particle 
case are magnified by a factor of 10 and 100, respectively. 

at THz waveband is larger than that in optical frequency. 
The spot size of the laser beams with appropriate inci- 
dent angles is taken as l x x l y . Thus the power of the 
THz signal propagating along the in-plane x-direction is 
roughly 

P T H Z «w 2 |x diff | 2 ^ J D(8/ !/ c 3 eoe 1 / 2 n 2 ) * P^, (5) 

where Pi and Pi are the power of the two input laser 
beams, N is the number of the superlattice periods, and 
c is the vacuum light velocity. To be specific, we use 
following realistic parameters: N = 50, P\ = Pi = 0.1 
W, and l y = 1 mm. 

In Fig. |21 the difference susceptibility and the THz 
emission power are plotted against the output frequency 
both for a fixed static field and for doubly resonant condi- 
tion. For comparison, the results without Coulomb inter- 
action are also shown for the fixed field strength. The ex- 
citonic effect enhances the emission power by more than 
two orders of magnitude. Under doubly resonant con- 
dition, the power of THz emission is several /iW, and 
the THz electric field strength is larger than 1 kV/m. 
The efficiency of converting near infrared photons to THz 
photons is around 1%. Because the doubly resonant con- 
dition can always be realized by simultaneously adjusting 
the bias voltage and laser frequencies, the THz emission 
power is not sensitive to the frequency difference in the 
range of 0.5-6 THz, which, under fixed static field, would 



otherwise decreases rapidly as the frequency difference 
goes away from the Bloch oscillation frequency. 

Negligible absorption of THz signals is another advan- 
tage of using biased superlattices for difference-frequency 
over other doubly-resonant schemes, such as that in 
doped quantum wells. For parameters in the example 
above, the optically excited carrier density, with contin- 
uous wave operation mode and 1 ns _1 recombination rate 
assumed, is of the order of 10 9 cm~ 2 per quantum well, 
and there is basically no free carriers outside the laser 
spot, so the THz signal can propagate without significant 
absorption by electrons. Neither is the signal absorption 
by optical phonons crucial in the frequency range consid- 
ered. 

In conclusion, the in-situ tunable THz emission from 
difference-frequency in biased superlattices is quite ef- 
ficient owing to the tunable doubly-resonant condition, 
the excitonic enhancement, and negligible absorption of 
signals. 
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